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Abstract 

Homodyne tomography — i. e. homodyning while scanning the local oscillator 
phase — is now a well assessed method for "measuring" the quantum state. In this 
paper I will show how it can be used as a kind of universal detection, for mea- 
suring genèric field operators, however at expense of some additional noise. The 
general class of field operators that can be measured in this way is presented, and 
includes also operators that are inaccessible to heterodyne detection. The noise 
from tomographical homodyning is compared to that from heterodyning, for those 
operators that can be measured in both ways. It turns out that for some operators 
homodyning is better than heterodyning when the mean photon number is suffi- 
ciently small. Finally, the robustness of the method to additive phase-insensitive 
noise is analyzed. It is shown that just half photon of thermal noise would spoil the 
measurement completely. 



1 Introduction 



Homodyne tomography is the only viable method currently known for determining the detailed 
state of a quantum harmònic oscillator — a mode of the electromagnètic field. The state mea- 
surement is achieved by repeating many homodyne measurements at different phases (fi with 
respect to the local oscillator (LO). The experimental work of the group in Eugene-Oregon Q 
undoubtedly established the feasibility of the method, even though the earlier data analysis were 
based on a filtered procedure that affected the results with systematic errors. Later, the theoret- 
ical group in Pavia-Italy presented an exact reconstruction algorithm Q, which is the method 
currently adopted in actual experiments (see, for example, Refs. || and The reconstruction 
algorithm of Ref. || was later greatly simplified 0, so that it was possible also to recognize the 
feasibility of the method even for nonideal quantum efficiency ry < 1 at the homodyne detector, 
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and, at the same time, establishing lower bounds for rj for any given matrix representation. After 
these first results, further theoretical progress has been made, understanding the mechanisms 
that underly the generation of statistical errors ||, thus limiting the sensitivity of the method. 
More recently, for 77 = 1 non trivial factorization formulas have been recognized [[/], ^] for the 
"pattern functions" j|] that are necessary to reconstruct the photon statistics. 

In this paper I will show how homodyne tomography can also be used as a method for 
measuring genèric field operators. In fact, due to statistical errors, the measured matrix elements 
cannot be used to obtain expectations of field operators, and a different algorithm for analyzing 
homodyne data is needed suited to the particular field operator whose expectation one wants to 
estimate. Here, I will present an algorithm vàlid for any operator that admits a normal ordered 
expansion, giving the general class of operators that can be measured in this way, also as a 
function of the quantum efficiency r/. Hence, from the same bunch of homodyne experimental 
data, now one can obtain not only the density matrix of the state, but also the expectation 
value of various field operators, including some operators that are inaccessible to heterodyne 
detection. However, the price to pay for such detection flexibility is that all measured quantities 
will be affected by noise. But, if one compares this noise with that from heterodyning (for those 
operators that can be measured in both ways), it turns out that for some operators homodyning 
is less noisy than heterodyning, at least for small mean photon numbers. 

Finally, I will show that the method of homodyne tomography is quite robust to sources of 
additive noise. Focusing attention on the most common situation in which the noise is Gaussian 
and independent on the LO phase, I will show that this kind of noise produces the same effect of 
nonunit quantum efficiency at detectors. Gener alizing the result of Ref. ||], I will give bounds 
for the overall rms noise level below which the tomographical reconstruction is still possible. 
I will show that the smearing effect of half photon of thermal noise in average is sufficient to 
completely spoil the measurement, making the experimental errors growing up unbounded. 



2 Short up-to-date review on homodyne tomography 

The homodyne tomography method is designed to obtain a general matrix element {ip\ g\(p) in 
form of expectation of a function of the homodyne outcomes at different phases with respect to 
the LO. In equations, one has 

f w d<p f +00 

(ip\g\(p) = — / dxp{x;4>) /^(x; ç!>) , (1) 

JO 7T J-<x> 

where p(x; 4>) denotes the probability distribution of the outcome x of the quadrature x^ = 
i (a^e 1 ^ + ae'^j of the field mode with particle operators a and a) at phase (j) with respect 
to the LO. Notice that it is sufficient to average only over cj) 6 [0, 7r], due to the symmetry 
x<^ +7r = —Xef,. One wants the function /^(x; (fi) bounded for all x, whence every moment will be 
bounded for any possible (a priori unknown) probability distribution p(x;<j)). Then, according 
to the central-limit theorem, one is guaranteed that the integral in Eq. ([]]) can be sampled 
statistically over a sufficiently large set of data, and the average vàlues for different experiments 
will be Gaussian distributed, allowing estimation of confidence intervals. If, on the other hand, 
the kernel f^{x\ 4>) turns out to be unbounded, then we will say that the matrix element cannot 
be measured by homodyne tomography. 

The easiest way to obtain the integral kernel f^[x\ <f>) is starting from the operator identity 

q= í _Tr(èe- SaW )e- aat+Sa (2) 
J vr 
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which, by changing to polar variables a = (i/2)ke t ^, becomes 

Tr(ge ikx ^) e~ lkx ^ . (3) 



n d<l> f +0 ° dk\k\ itó 



vr 

Equation (g) is nothing but the operator form of the Fourier-transform relation between Wigner 
function and characteristic function: it can also be considered as an operator form of the Moyal 
identity 

(l 2 7 

— {k\D\z)\m){l\D{z)\n) = (k\n)(l\m) . (4) 

The trace-average in Eq. @ can be evaluated in terms of p(x, <f>), using the complete set {|;r}</>} 
of eigenvectors of x^, and exchanging the integrals over x and k. One obtains 

dxp{x\ (j))K{x — xq) , (5) 
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where the integral kernel K{x) is given by 

K(x) = --p4t = - lim -Re- , (6) 

w 2 x 2 e->o+ 2 (x + ie) 2 ' v ; 

P denoting the Cauchy principal value. Taking matrix elements of both sides of Eq. (||) between 
vectors ip and (p, we obtain the sampling formula we were looking for, namely 

fT r+co 

(ip\g\<p) = — / dxp(x;4>)(ip\K(x-x <j) )\(p). (7) 
Jo vr J_oo 

Hence, the matrix element {^\g\ip) is obtained by averaging the function ff v (x; ep) = (ip\K(x — 
xj))\(p) over homodyne data at different phases ep. As we will see soon, despite K{x) is unbounded, 
for particular vectors tp and 92 in the Hilbert space the matrix element {ip\K(x — Xé)\tp) is 
bounded, and thus the integral ([?]) can be sampled experimentally. 

Before analyzing specific matrix representations, I recali how the sampling formula (0) can 
be generalized to the case of nonunit quantum efhciency. Low efficiency homodyne detection 
simply produces a probability p v (x;(p) that is a Gaussian convolution of the ideal probability 
p(x;cp) for 7] = 1 (see, for example, Ref. fl~0|| ). In terms of the generating functions of the 
x^-moments one has 

+00 / 1 „ \ r+oo 

dxp v (x;<P)e ikx = exp — -^fc 2 ) / dxp(x; 0)e ikx . (8) 

-00 V OT] J J — 00 

Upon substituting Eq. ^ into Eq. (|3|) , and by following the same lines that lead us to Eq. (|||) , 
one obtains the operator identity 

f n dep f +oc 

g= — / dxp v (x;(j))K rt (x - , (9) 



vr 



where now the kernel reads 



KJx) = -Re / + dkk exp (- — -k 2 + ikx] . (10) 
2 Jo \ 8r] ) 

The desired sampling formula for {ip\g\(p) is obtained again as in Eq. (|7|), by taking matrix 
elements of both sides of Eq. (|Ï0|), Notice that now the kernel K v (x) is not even a tempered 
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distribution: however, as we will see immediately, the matrix elements of K„(x—x,p) are bounded 
for some representations, depending on the value of rj. The matrix elements (iplK^x — à^)\ip) 
are bounded if the following inequality is satisfied for all phases <p £ [0, ir] 



where e 2 (ó) is the harmònic mean 



2 11 

^W)~4Ï^) + W T: ' (12) 



and £y(((>) is the "resolution" of the vector \v) in the x^-representation, namely: 

„2 



(x|f)| 2 ~ exp 



x 



24(0) 



(13) 



In Eq. (13) the symbol ~ stands for the leading term as a function of x, and {x)^ = e iaa ^\x) 
denote eigen-ket of the quadrature x^ for eigenvalue x. Upon maximizing Eq. (lïï|) with respect 
to <j) one obtains the bound 

^ > 7T772 > e2 = J^M^)} ■ ( 14 ) 

1 + 4:E Z </>e[0,7r] 

One can easily see that the bound is r/ > 1/2 for both number-state and coherent-state repre- 
sentations, whereas it is rj > (1 + s 2 ) -1 > 1/2 for squeezed-state representations with minimum 
squeezing factor s < 1. On the other hand, for the quadrature representation one has r\ > 1, 
which means that this matrix representation cannot be measured. The value rj = 1/2 is actually 
an absolute bound for all representations satisfying the "Heisenberg relation" e(<p)e((j) + § ) > \ 
with the equal sign, which include all known representations (for a discussion on the existence 
of exòtic representations sec Rcf. [11]). Here, I want to emphasize that the existence of such a 
lower bound for quantum efficiency is actually of fundamental relevance, as it prevents measur- 
ing the wave function of a single system using schemes of weak repeated indirect measurements 
on the same system [12|. 

At the end of this section, from Ref. || I report for completeness the kernel (n\K(x — x$)\rri) 
for matrix elements between number eigenstates. One has 



{n\KJx - x^)\n + d) = e -^ 2 ^+2 / w " e -" 2 * 2 (15) 

y (n + d)\ 

j2tLí n + d \ (2„ + d + l)!K 2 "Re U-i) d D_ [2v+d+2) (-2ÍKx)} , 

— LS . \ I 1/ LS / ^ 



u=0 



where k = \/r//(2r] — 1), and D a (z) denotes the parabòlic cylinder function. For rj = 1 the 
kernel factorizes as follows j/], |8| 

(n\K(x — x$)\n + d) 



= e 1 [2xu n (x)u n+(i (2;) - y/n + lu n+ i(x)v n+d (x) - \[m + ïu n (x)v n+d +i(x)} , (16) 

where u n (x) and v n (x) are the regular and irregular energy eigen-functions of the harmònic 
oscillator 
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3 Measuring genèric field operators 



Homodyne tomography provides the maximum achievable information on the quantum state, 
and, in principle, the knowledge of the density matrix should allow one to calculate the expec- 
tation value (O) = Tr[Og] of any observable O. However, this is generally true only when one 
has an analytic knowledge of the density matrix, but it is not true when the matrix has been 
obtained experimentally. In fact, the Hilbert space is actually infinite dimensional, whereas 
experiment ally one can achieve only a finite matrix, each element being affected by an experi- 
mental error. Notice that, even though the method allows one to extract any matrix element 
in the Hilbert space from the same bunch of experimental data, however, it is the way in which 
errors converge in the Hilbert space that determines the actual possibility of estimating the trace 
Tr[O^J. To make things more concrete, let us fix the case of the number representation, and 
suppose we want to estimate the average photon number (a)a). In Ref. |Ï3| it has been shown 
that for nonunit quantum efficiency the statistical error for the diagonal matrix element (n\g\n) 
diverges faster than exponentially versus n, whereas for w = 1 the error saturates for large n 
to the universal value e n = y/2/N that depends only on the number N of experimental data, 
but is independent on both n and on the quantum state. Even for the unrealistic case r\ = 1, 
one can see immediately that the estimated expectation value (a^a) = J2n=o n Qnn based on the 
measured matrix elements g nn , is not guaranteed to converge versus the truncated-space dimen- 
sion H, because the error on g nn is nonvanishing versus n. Clearly in this way I am not proving 
that the expectation (a^a) is unobtainable from homodyne data, because matrix errors con- 
vergence depends on the chosen representation basis, whence the ineffectiveness of the method 
may rely in the data processing, more than in the actual information contained in the bunch of 
experimental data. Therefore, the question is: is it possible to estimate a genèric expectation 
value (O) directly from homodyne data, without using the measured density matrix? As we will 
see soon, the answer is positive in most cases of interest, and the procedure for estimating the 
expectation (O) will be referred to as homodyníng the observable O. 

By homodyníng the observable O I mean averaging an appropriate kernel function 7£.[Ó](ic; (f)) 
(independent on the state g) over the experimental homodyne data, achieving in this way the 
expectation value of the observable (O) for every state g. Hence, the kernel function 7?.[0](a;; 4>) 
is defined through the identity 

(ó) = r# í + °° dxp{x;tf>)n[Ó](x;<l>). (18) 

JO 7T J-oo 

From the definition of 1Z[Ò](x; <p) in Eq. flï~8|), and from Eqs. (Q) and (Q) — which generally hold 
true for any Hilbert-Schmidt operator in place of g — one obtains 

6= f °° dx <ji[Ò]{x;iï\x) M (x\ , (19) 



lo vr 

with the kernel lZ[0](x; (j)) given by 

K{Ò](x; <P) = Tr[ÒK(x - X(f> )] , (20) 
and K(x) given in Eq. (0). The validity of Eq. (pÜ[), however, is limited only to the case of a 



Hilbert-Schmidt operator O, otherwise it is ill defined. Nevertheless, one can obtain the explicit 
form of the kernel 1Z[0](x; (fi) in a different way. Starting from the identity involving trilinear 



products of Hermite polynomials [14| 



m+n+fe 1 

+°° _ 2 2 2 7r2/c!m!n! 

dx e x H k (x) H m (x) H n (x) = — — , for k + m + n = 2s even , (21) 

(s — k)\(s — m)\(s — n)\ 



-oo 
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Richter proved the following nontrivial formula for the expectation value of the normally ordered 
field operators [15] 



{a) n a m ) 

lO 7T 

which corresponds to the kernel 



+00 u 

i(m-n)<j> n n+m 



dxp(x; 4>)e l 



{V22 



v /2«+ m ("-+" 1 ) ' 



K[a) n a m ](x-(t)) 



J(m-n)(j> Hn+ m 



(V2x) 



(22) 



(23) 



This result can be easily extended to the case of nonunit quantum efficiency 77 < 1, as the 
normally ordered expectation (a^ n a m ) just gets an extra factor 7j2^ n+m \ Therefore, one has 



K v [aï n a m ](x;4>) =é (m ~ n) 



(24) 



V(2í?) n +m( n+ m) > 

where the kernel TZ v [0](x; (f>) is defined as in Eq. (|Ï8|), but with the experimental probability 
distribution p v (x;4>). From Eq. (p4| ) by linearity on can obtain the kernel TZrj[f](x; 4>) for any 
operator function / that has normal ordered expansion 



/ = /(a,at)= /, 



( n ) r A n n m 



nm=0 



From Eq. fl24| ) one obtains 
" H s (V2x) 



n, n [f]{x-^) 

where 



E 

s=0 



í!(2? 7 ) s / 2 



E /, 

nm=0 



(n) i(m-n) 



)l 



!m!<5- 



n+m,s 



E 

s=0 



H s {V2x)i s d s 



(25) 



J·[/]( U ;0),(26) 

d=0 



nm=ü 

Continuing from Eq. (p6[) one obtains 



_l v \n+m e i{m-n) 



_ , i, , . . ( 1 d? 2ix d " 

/ = exp — —2 + — — 

V 2?7 mr ai> 



u=0 



and finally 



+00 



ífu> 



=« -i w V[/](«; + 2^/^;. 



(27) 



(28) 



(29) 



Hence one concludes that the operator / can be measured by homodyne tomography if the 
function J r [/](í;; </>) in Eq. (27) grows slower than exp(— r]v 2 /2) for v — > 00, and the integral 
in Eq. fl29| ) grows at most exponentially for x — > 00 (assuming p(x; <f>) goes to zero faster than 
exponentially at x — > 00). 

In Table [ï] I report the kernel 7^. ?? [Ó](x;0) for some operators Ó. One can see that for 
the raising operator è+ the kernel diverges at 77 = l/2 + , namely it can be measured only for 
ï] > 1/2. The operator W s in the same table gives the generalized Wigner function W s (a, a) for 
ordering parameter s through the identity W s (a,a) = Tr[D(a)gD^ (a)W s ]. From the expression 
of TZ v [W s ](x; <p) it follows that by homodyning with quantum efficiency r] one can measure the 
generalized Wigner function only for s < 1 — r/ _1 : in particular, as already noticed in Refs. ||, 
the usual Wigner function for s = cannot be measured for any quantum efficiency [in fact one 
would have Hi[D\a)W D{a)){x; (/)) = K[x - Re(ae~^)], with K(x) unbounded as given in Eq. 
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o 



7^[0](^) 



a tn a m 



J(m-n)<j> Hn+m 



2x) 



2e l(p x 



3 2 ^(4x 2 - 1) 



(3) 



i — r 



(4) 



2x 



(aJaY 



|x 4 - 2x 2 



(5) 



exp[-^(ae^) 2 + ^ae l 
1 + ge~ 2i * ~ 



exp[-^(ae-^) 2 - ^ae" 1 * 



(6) 



Dt( 



ei 



7.T ri/(T 



+ 



1 + 

OL 



(7) 



e + = a< 



\/ï+ota 



2xe 



v 7 ^ 



+oo 



ric 



r$ 2, 



.1 



2' 1 + 2- 



2jj 



(8) 



W - 2 ( s +l 



di- 



2e" 



7r(l — s) 



cos 2, 



2í 



(1- S ) 



(9) 



|n + d)(n| 



(n|.K"(x — X0)|n + d) in Eqs. (|Ï5 ) and (1 



Table 1: Kernel TZ v [Ò](x- 1 cj>), as defined in Eq. (pt), for some operators Ò. [The symbol <i>(a,ò;a;) 
denotes the customary confluent hypergeometric function.] 



3.1 Comparison between homodyne tomography and hetero- 
dyning 

We have seen that from the same bunch of homodyne tomography data, not only one can recover 
the density matrix of the field, but also one can measure any field observable / = /(a, cd) 
having normal ordered expansion / = f^ n \a,a)) = Y^ m =o ímn(^ n a m and bounded integral in 
Eq. (pG|) — this holds true in particular for any polynomial function of the annihilation and 
creation operators. This situation can be compared with the case of heterodyne detection, 
where again one measures general field observables, but admitting anti-normal ordered expansion 
/ = f^ a H a i °0 = Y^ m =o fnma m a^ n , in which case the expectation value is obtained through the 
heterodyne average 

</> = / — f {a) (a,a)(a\è\a) . (30) 
J vr 

For j] = 1 the heterodyne probability is just the Q-function Q(a,a) = ^(a\g\a), whereas for r) = 
1 it will be Gaussian convoluted. As shown by Baltin |Ï6|| , generally the anti-normal expansion 
either is not defined, or is not consistent on the Fock basis, namely f( a \a, a))\n) has infinite 
norm or is different from f(a,a))\n) for some n > 0. In particular, let us focus attention on 
functions of the number operator f{a)a) = Y^Ï^=o c l( a ^ a ) 1 , f^(a^a) = Z^o c | c^ l a l , f^(a)a) = 
T^o^a 1 ^ 1 . Baltin has shown that |Ï6j 



1 = nL dX9(x){e _1) = £ m- k y. ' 

^^(A^^^ /lV 1 . (31) 



+0 ° dx s iXx 



5(A) = / -f(x)e 
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From Eqs. ( |3l| ) one can see that the normal ordered expansion is always well defined, whereas 
the anti-normal ordering needs extending the domain of / to negative integers. However, even 
though the anti-normal expansion is defined, this does not mean that the expectation of /(a^a) 
can be obtained through heterodyning, because the integral in Eq. ( |30| ) may not exist. Actually, 
this is the case when the anti-normal expansion is not consistent on the Fock basis. In fact, for 
the exponential function f(a^a) = exp(— fiaJa) one has f( a \\a\ 2 ) = e ít exp[(l — e^)|a| 2 ]; on the 
Fock basis f^ a \a) a)\n) is a binomial expansion with finite convergence radius, and this gives the 
consistency condition |1— é í \ < 1. However, one can take the analytic continuation corresponding 
for 1 — é 1 < 1, which coincides with the condition that the integral in Eq. (|30|) exists for any 
state g (the Q-function vanishes as exp(— |q| 2 ) for a — > oo, at least for states with limited photon 
number). This argument can be extended by Fourier transform to more general functions f(a'a), 
leading to the conclusion that there are field operators that cannot be heterodyne-measured, 
even though they have well defined anti-normal expansion, but the expansion is not consistent on 
the Fock basis. As two examples, I consider the field operators è+ and W s in Table p]. According 
to Eqs. (|3Ï|) it follows that the operator è+ does not admit an anti-normal expansion, whence 
it cannot be heterodyne detected. This is in agreement with the fact that according to Table 
|H we can homodyne è+ only for r\ > 1/2, and heterodyning is equivalent to homodyning with 
effective quantum efficiency rj = 1/2 (which corresponds to the 3 dB noise due to the joint 
measurement |Ï7]]). The case of the operator W s is different. It admits both normal-ordered and 

anti-normal-ordered forms: W s = n ^i_ s ^ '■ exp jr^o^o) := ~ 7r (i+ s ) : ex P (ï^°^ a ) '-A, where 
: . . . : denotes normal ordering and : . . . \a anti-normal. However, the consistency condition for 
anti-normal ordering is 2/(s + 1) < 1, with s < 1, which implies that one can heterodyne W s for 
s > —1, again in agreement with the value of s achievable by homodyne tomography at r\ = 1/2. 

Now I briefly analyze the additional noise from homodyning field operators, and compare 
them with the heterodyne noise. For a complex random variable z = u + iv the noise is given 
by the eigenvalues = \z\ 2 — \z\ 2 db \z 2 — ~z 2 \ of the covariance matrix. When homodyning 

the field, the random variable is z = 2e %< ^x Jïg] and the average over-line denotes the doublc 
integral over x and <j) in Eq. (18). From Table (||) one has ~z = (a), z 2 = (a 2 ), \z\ 2 = 2{a)a) + 1, 



e 2 ^ = 0§. In this way one finds that the noise from homodyning the field is A^]ja] 
1 + 2(a)a) — \(a)\ 2 ± |(« 2 ) — ( a ) 2 |- On the other hand, when heterodyning, z becomes the 
heterodyne output photocurrent, whence z = (a), z 2 = (a 2 ), \z\ 2 = (a'a) + 1, and one has 
Nf^[a] = 1 + (a)a) — \ (a)\ 2 ± |(a 2 ) — (a) 2 \, so that the tomographical noise is larger than the 
heterodyne noise by a term equal to the average photon number, i. e. 

Njg l [o]=Jv£ ) [o] + (oto>. (32) 

Therefore, homodyning the field is always more noisy than heterodyning it. On the other 
hand, for other field observables it may happen that homodyne tomography is less noisy than 
heterodyne detection. For example, one can easily evaluate the noise A^ om [n] when homodyning 
the photon number n = a^a. The random variable corresponding to the photon number is 
v(z) = |(|.z| 2 — 1) = 2x 2 — |, and from Table |ï| we see that the noise Nh ortl [n\ = Au 2 (z) 
can be written as iV/j 0r7l [n] = (An 2 ) + ^(n 2 + n + 1) ]Ï3] . When heterodyning the field, the 
random variable corresponding to the photon number is v(z) = |z| 2 — 1, and from the relation 
|z| 4 = {a) 2 a 2 ) one obtains Nh e t[n] = Ais 2 (z) = (An 2 ) + (n + 1), namely 

N hom [n] = N het [n] + ^(n 2 - n - 1) . (33) 

We thus conclude that homodyning the photon number is less noisy than heterodyning it for 
sufficiently low mean photon number (n) < |(1 + 
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4 Homodyne tomography in presence of additive phase- 
insensitive noise 



In this section I consider the case of additive Gaussian noise, in the typical situation in which 
the noise is phase-insensitive. This kind of noise is described by a density matrix evolved by the 
master equation 



dtg(t) = 2 AL[a)]+BL[a] g(t) 



(34) 



where L[c] denotes the Lindblad super-operator L[c]g = cgc* — ^{&c, g]+. Due to the phase 
invariance L[ae~ l< ^] = L[a] the dynamical evolution does not depend on the phase, and the 
noise is phase insensitive. From the evolution of the averaged field (a) ou t = Tr[a£(í)] = g{a)i n = 
Tr[a^(0)] with g = exp[(A — B)t], we can see that for A > B Eq. (|34| ) describes phase-insensitive 
amplification with field-gain g, whereas for B > A it describes phase-insensitive attenuation, 
with g < 1. Concretely, for A > B Eq. (|34| ) models unsaturated parametric amplification 
with thermal idler [average photon number fh = B/(A — B)], or unsaturated làser action [A 
and B proportional to atòmic populations on the upper and lower lasing levels respectively]. 
For B > A, on the other hand, the same equation describes a field mode damped toward 
the thermal distribution [inverse photon lifetime T = 2{B — A), equilibrium photon number 
fh = A/(B — A)], or a loss g < 1 along an optical fiber or at a beam-splitter, or even due to 
frequency conversionpOj]. The borderline case A = B leaves the average field invariant, but 
introduces noise that changes the average photon number as {a)a) out = {a)a)i n + n, where 
n = 2At. In this case the solution of Eq. (S3) can be cast into the simple form 



g(t) 



d 2 (3 



nn 



exp 



(-\(3\ 2 /n)D((3)Q(0)D\(3). 



(35) 



This is the Gaussian displacement noise studied in Refs-plj |22| and commonly referred to as 
"thermal noise" [regarding the misuse of this terminology, see Ref. [2^| ], which can be used to 
model many kinds of undesired environmental effects, typically due to linear interactions with 
random classical fluctuating fields. 

Eq. ( |34|) has the following simple Fokker-Planck differential representation |^] in terms of 
the generalized Wigner function W s (a, a) for ordering parameter s 



d t W s {a,a;t) = Q(d a a + d & a) + 2D s d 2 a a W s (a,a;t) , 



(36) 



where Q = B — A and 2D S = A + B + s(A — B). For nonunit quantum efficiency r\ and after a 
noise-diffusion time í the homodyne probability distribution p J1 (x;cj);t) can be evaluated as the 
marginal distribution of the Wigner function for ordering parameter s = 1 — namely 



/+ oo , . 

dyW 1 _ r) -x Ux + iy)e t *, (x - iy)e _, *;t) 
-oo v ' 

The solution of Eq. ([36]) is the Gaussian convolution [ 23 



W s (a,a;t) 



[d 2 /3 


\a-g(3\ 2 ' 


/^f eXP 


[ s 2 



W a (J3,P;0) , 



5 2 = ^(l-e~ 2 ^) 



and using Eq. (|37| ) one obtains the homodyne probability distribution 

dx' 



p v (x;(j);t) = e Qt 



-oo j2irA 2 1 _ rj _ 1 



exp 



(x' - g 1 x) 2 
i— i) 1 



p v {x';4>) . 



(37) 



(38) 



(39) 
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where = ^d^^i-n- 1 ' It is easy to see that the generating function of the í^-moments with 
the experimental probability p„(x;(fi;t) can be written in term of the probability distribution 
p(x; (fi) for perfect homodyning as follows 

/+oo /i 1 \ r+oo 

^ dx Pv (x; (fi; t)e lkx = exp [~-g 2 ^k 2 - -g-Wj / ^ dxp(x; (fi)e^ kx , (40) 

Eq. ( |40| ) has the same form of Eq. (||) , but with the Fourier variable k multiplied by g and with 
an overall effective quantum efficiency rj* given by 

C 1 = rT 1 + 4A^ = 5 - V 1 + -^jid- 2 ~ 1) • (41) 

On the other hand, following the same lines that lead us to Eq. (H), we obtain the operator 
identity 

f n dó f +oc 

q = q(Q)= — / dxp v ^(x;(t>;t)K~(g~x-X6), (42) 

JO 7T J-oo 

which also means that when homodyning the operator O one should use lZ rit (g~ 1 x](fi) in place 
of Hn(x; (fi), namely, more gener ally, one needs to re-scale the homodyne outcomes by the gain 
and use the effective quantum efficiency 77* in Eq. (41). In terms of the gain g and of the 
input-output photon numbers, the effective quantum efficiency reads 

V* 1 = rT X + g- 2 {2{a ] a) out + O - (2(a f a) m + r?" 1 ) . (43) 

In the case of pure displacement Gaussian noise (A = B), Eq. (|43|) becomes 

r]- 1 = r]- 1 + 2n , (44) 

which means that the bound r/* > 1/2 is surpassed already for n > 1: in other worlds, it is just 
sufficient to have half photon of thermal noise to completely spoil the tomographic reconstruc- 
tion. 
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